FINAL-STATE RADIATION AND LINE-SHAPE DISTORTION 
IN RESONANCE PAIR PRODUCTION 

W. Beenakker*\ F.A. Berends and A. P. Chapovsky^^ 

Instituut-Lorentz, University of Leiden, The Netherlands 



ABSTRACT 



In this letter it is shown how final-state QED corrections to the production of a pair of reso- 
nances can distort the line shape of such a resonance in a sizeable way. This effect depends 
on the definition of the line shape and can reach up to 30%, depending on the final state. 
The mechanism is first displayed for a particular case of ZZ production, for which an exact 
and approximate treatment can be given. The approximate method is then applied to ly-pair 
production. In addition some simple rules of thumb are given for accurately estimating the 
characteristic distortion effects, like the mass shift and peak reduction. 
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1 Introduction 



As is well known the Z line shape as measured in e"'"e^ ^ Z — 2/ is distorted due to initial- 
state radiation (ISR). Without ISR the total cross-section a{s) as a function of the square of the 
centre-of-mass energy s gives the line shape. With ISR the centre-of-mass energy available to 
produce a Z boson changes and, as a consequence, so does the shape of the total cross-section 
cr(s). Experimentally the latter is measured. If one would measure the square of the modified 
centre-of-mass energy s', one would determine cr(s') and thereby the pure line-shape. It should 
be noted that final-state radiation (FSR) only marginally corrects the overall size of cr(s), but 
not its shape. Therefore FSR is less relevant for the usual Z line-shape measurement. 

When one produces two resonances, or one resonance and a stable particle, the line shape of 
such a resonance will be measured from the invariant-mass distribution of its decay products. 
Examples are pair production of W bosons, Z bosons, tt or HZ . Depending on how one 
measures the invariant-mass distribution of the decay products of the particular resonance, one 
finds the pure line shape or a distorted one. This time also FSR can cause the distortion. 

It is the main purpose of this letter to point out that such a FSR-induced distortion can 
arise. For exhibiting the effect we take an example for which we can perform both exact and 
approximate calculations. An ideal example is the double-resonance process 

v^v^ ^ ZZ ^ e+e^z/^z/^. (1) 

Here QED corrections apply only to the decay Z e+e~ and not to the other Z decay or the 
initial state. When the Z line shape is obtained from measuring the invariant-mass distribution 
of the e^e~ pair, FSR will distort it in a way reminiscent of the usual ISR distortion in single 
Z production. The virtue of the example is threefold. In the first place, process (|I]) is free 
of the gauge-invariance problems that are inherent in the production of unstable particles. 
This holds in spite of the fact that we have left out all non-double-resonant mechanisms for 
producing the e^e'v^-VT final state. Secondly, the QED radiative corrections only lead to FSR. 
So, the effect of FSR on the line shape can be studied without the additional presence of ISR 
phenomena. Thirdly, the effect can be calculated exactly. In more realistic examples, involving 
for instance e^e~ initial states, Z-pair production with both Z bosons decaying into charged 
particles, or ly-pair production, additional classes of QED radiative corrections emerge, like 
ISR [0] or non-factorizable interference corrections 0. Moreover, in order to avoid gauge- 
invariance problems, the QED corrections often have to be calculated in an approximation, 
which for instance restricts the calculation to the leading logarithmic corrections and/or the 
leading terms in a pole-scheme expansion around the resonances 0, ^. Nevertheless the FSR 
distortion of the line shape will still be one of the main features. In these more complicated 
cases Monte Carlo studies including radiative corrections would be needed. Here we focus 
exclusively on the line-shape deformation and its impact on the determination of the resonance 
mass. 

Although we start with reaction (|1]), we shall also comment on the more realistic case of 
the W line shape at LEP2. 
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2 The Z^-pair example: exact calculation 



For process (|l|) we first consider the Born approximation, to which two double-resonant dia- 
grams contribute. After integration over the Z production angle and the fermion decay angles, 
one obtains 



daoiMf,Ml 



A2(M|) 



|Di(M2)|2 \D,{MI)\^ 



F{M^,M. 



2>i 



(2) 



where and M| denote the invariant masses of the e^e and Vt-Vj. pairs, respectively. The 
spin-averaged production cross-section takes the form 
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with A the Kallen function 

A = + + - 2{sMl + sMl + M^M^), 
which is in agreement with the literature [H. The decay parts are given by 
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whereas the resonance shapes are dominated by 
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Di,2(Mj2)=M22-M| + ^Mj2 
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Note that we have used the standard LEPl representation in the above formulae, involving Gp 
and the effective couplings of the Z boson to leptons {gyi-, Qai) and neutrinos (g^). 



Applying virtual and soft photonic corrections to (0) yields 
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Here we have defined the soft photons in the rest frame of the Z: < eMi/2 -C Tz- Photon 
bremsstrahlung involving more energetic photons introduces an explicit dependence on the 
photon energy E^, resulting in a distribution in the invariant masses of both the e^e~ pair 
{Mf) and the e+e~7 system [Mf = virtuality of the Z boson): 
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where 
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When correction (|^) is combined with (P) and an integration over Mf is performed, the cor- 
rection to da^/ [dMl dM2) takes on the form of the usual FSR factor 1 + 3a/(47r). This is 
in agreement with the KLN theorem, which imphes that the large logarithmic contributions 
(oc L, L') vanish upon summation (integration) over all degenerate final states. So, the reso- 
nance shape is not deformed when one measures the Ml distribution, i.e. the invariant-mass 
distribution of the e"'"e~7 system. 

In our special example this choice of distribution is, of course, the natural one. However, in 
more realistic processes it is in general unclear whether the photon is radiated from the initial 
state, the unstable particles, or the final state. This introduces the freedom to either choose 
Ml or Ml for the definition of the invariant-mass distribution of the unstable particle. 

If one measures the Ml distribution one will find a distorted line shape. The reason is 
that (ID now has to be integrated over values ranging from Ml to — M2Y ■ This causes 
the Ml line shape to receive contributions from effectively higher Z-boson virtualities. This is 
to be compared with the single- Z-production case where the ISR-corrected line shape receives 
contributions from effectively lower Z-boson virtualities. Due to the fact that roughly speaking 
the resonance shape is symmetric around the resonance mass, one expects now a distortion of 
the resonance shape that is approximately the LEPl distortion reflected with respect to the 
resonance mass. 

The bremsstrahlung contribution to the line shape da / [dMl dM2) arising from (^ reads 
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where the function F is defined in (0) and Cmax = {^/s — M2)^/Mf . Combining (|^) and (|TT]) 
gives the 0{a) correction to the line shape. Just like at LEPl it will be necessary to resum 
the soft corrections, as will become clear from the discussion in the following section. Based on 
LEPl experience |]l| a suitable expression for this resummation is given by 
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with 
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Figure 1: The FSR-induced distortion of the hne shape da/i^dMldM^) corresponding to pro- 
cess (0) for M2 = Mz- Centre-of-mass energy: -^i = 200 GeV. 

In Fig. 0we display the FSR-induced distortion effects on the hne shape da/{dMf c/M|) for 
a centre-of-mass energy of a/s = 200 GeV and a fixed invariant mass M2 = Mz- However, the 
actual distortion phenomena do not depend on the precise value of M2. The parameter input 
used in the numerical evaluation is: 

Mz = 91.1867 GeV, Tz = 2.4948 GeV, Gp = 1.16639 x 10"^ GeV"^ = 0.51099906 MeV, 
= 137.0359895, = 0.50125, gve = -0.03681, gAi = -0.50112. 

The sizeable distortion effects are clearly visible, just as the importance of the soft-photon 
resummation. Compared with the Born line shape, the 0{a) (resummed) QED corrections 
induce a shift in the peak position of —199 MeV (—112 MeV) and a reduction of the peak 
height by 29% (26%). The size of these effects are a direct result of the non-cancellation of 
the leading logarithmic corrections, which can be understood from the observation that a fixed 
value for Mi makes it impossible to sum over all degenerate final states. Another noteworthy 
observation is the close similarity of the curves in Fig. |l] to the ones for the Z line shape at 
LEPl [|l|. As predicted, the two sets of curves are approximately related by reflection with 
respect to the Born peak position. 




3 The Z-paiv example: approximations 

As mentioned before, in order to calculate QED corrections to more realistic processes like W- 
pair production one in general has to resort to approximations. First of all, the fact that we are 
dealing with unstable (charged) particles introduces the problem of a gauge-invariant treatment 
of the finite- width effects . An appropriate way of handling this problem is by applying the 
pole scheme, i.e. by performing an expansion around the resonances. When it comes to 0{a) 
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corrections, it is sufficient to consider only the leading (double-pole) term in this expansion, 
leaving out terms that are formally suppressed by at least aLT / [ttM). From now on we 
will refer to this procedure as the double-pole approximation and indicate quantities that are 
calculated in this approximation by a bar. Note that the approximation only makes sense near 
the resonance of the unstable particle and sufficiently far above the production threshold of the 
underlying on-shell production process. The latter is caused by the direct relation between the 
double-pole residues and the on-shell production and decay processes. Secondly, as we have seen 
in the previous section, the leading logarithmic corrections (oc V) constitute the bulk of the 
FSR distortion effects. Moreover, these leading-log effects are universal and gauge invariant, 
being directly related to the collinear limit of photon radiation off light particles (like e^). 
In particular the universality property is appealing, since it implies that the description of 
the leading-log corrections does not depend on the specific features of the unstable particles 
and their photonic interactions. Therefore, it is worthwhile to further restrict the double- 
pole calculations to the leading logarithms. This additional approximation is referred to as 
the leading-log approximation (LLA). Before making any comments on processes like VT-pair 
production, we ffist concentrate on our Z-pair example and check the validity of the indicated 
approximations . 

We start off with the definition of the double-pole approximation. At Born level it amounts 

to 

oMM?^_ Ai(M|) A2(M|) 

dMl dMi ~ ^^^^^ |Di(M2)|2 \D,{Ml)r ^ ' 

with 

D,,2{Ml,) = Ml, - M| + tMzTz. (15) 

As mentioned before, the QED radiative corrections introduce an ambiguity in defining the 
invariant-mass distributions of the unstable particles. The most transparent way of illustrating 
this is by considering the case that the photon is radiated from such an unstable particle. For 
the description of the resonance before (after) radiation the natural choice of invariant mass 
involves the fermion pair with (without) the photon. The pole expansion can now in principle 
be performed around either resonance. In practice one has to choose one particular invariant 
mass for the distributions. For the purpose of studying FSR-induced distortion effects, we 
shall choose the e~^e~ invariant mass Mi in the following, although the e^e~7 invariant mass 
Ml would have been more natural for our special example (p. The corresponding double-pole 
approximation forces us to replace = (Mi by C^l; introducing an explicit dependence on 
the photon energy in the double-pole residues. This would even affect the (neutral) resonance- 
pair-production stage of the process. However, as can be verified explicitly, only semi-soft 
photons with energy = 0(Tz) ^ Mz contribute to the 0{a) corrected double-pole residues. 
As a result, ( can be effectively replaced by unity whenever possible, re-establishing the usual 
form of the double-pole approximation in terms of off-shell Breit-Wigner distributions and on- 
shell production/decay processes. The effects from hard photons {E^ ^ F^) are suppressed 
by at least Tz/Mz and are therefore neglected in the double-pole approximation. The picture 
underlying this phenomenon is that hard photons move the Z-boson virtuality (M^) far off 
resonance for near-resonance M^ values, resulting in a suppressed contribution to the line 
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shape. In fact, only the (soft) 1/(^ — 1) term in (pH]) contributes to the hne shape in the 
double-pole approximation. It should be noted that this very suppression of hard-photon effects 
serves as a postenon justification of the soft-photon resummation proposed in ([T2|) . 

With this observation in mind, the bremsstrahlung contribution ([TT|) takes the following 
form in double-pole approximation: 
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Here f3 can be derived from f3 by setting = M|. Note that the upper integration boundary 
Cmax has been extended to infinity, which is motivated by the fact that hard-photon effects 
are sufficiently suppressed. The remaining integral can be performed analytically. Combining 
with the virtual and soft corrections, which can be readily derived from (|^, we obtain in the 
double-pole LLA 
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In a way similar to the previous section the soft-photon corrections can be resummed, but this 
time the integral can be carried out explicitly in the double-pole LLA: 
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Having calculated the same quantities as in the previous section, we are now in the position 
to check the validity of the double-pole LLA. It turns out that the approximated results exhibit 
the same FSR distortions as the exact ones. Upon closer investigation, we observe for the 
0{a) (resummed) QED corrections a shift in the peak position of — 193MeV (— llSMeV) and 
a reduction of the peak height by 29% (26%). This is in excellent agreement with the distortion 
parameters of the exact calculation, proving the viability of the adopted approximations.^ 

Based on the results in the double-pole LLA, it is possible to derive simple and sufficiently 
accurate rules of thumb for the distortion parameters: 

• 0{a) corrections: the shift in the peak position AMf^'^'^ and the corresponding peak 
reduction factor ^p^'^'^ with respect to the Born line shape can be approximated by 

peak Vr/?rz/8 
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(19) 



^ For completeness we note that all curves in the double-pole LLA are displaced by a small amount with 
respect to the exact ones. This is caused by the fact that the Born results differ by subleading terms in the pole 
expansion. 
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• resummed corrections: now the distortion parameters read 

AMf"'^ ^ --l3Tzil + -) = -lllMeV, 
8 2 

This is in perfect agreement with the observed exact and double-pole distortion parameters. 
The analogy with the rules of thumb derived for the Z line shape at LEPl |^ confirms the 
relation between the FSR-induced distortion effects in double Z-resonance production and the 
ISR-induced distortion effects in single Z-resonance production at LEPl. 



4 Some comments on the W line shape at LEP2 

As has been shown in the previous section, the double-pole LLA constitutes a reliable framework 
for a gauge-invariant and universal description of FSR-induced distortion phenomena in double- 
resonance production. The essence of these phenomena is fully contained in the correction 



factor presented in (|18]), which apphes to each individual distorted Breit-Wigner distribution. 
For two distorted distributions the effect is hence multiplicative. Consequently, the reduction 
factor for a double-invariant-mass distribution is given by the product of the reduction factors 
for the individual single-invariant-mass distributions. However, the shift in the peak position 
does not change in the presence of more than one resonance; it only depends on the decay 
products of the unstable particle that is investigated. The only differences between process (|l]) 
and the more realistic process of VT-pair production at LEP2 are the resonance parameters 
{Mw = 80.22 GeV and Tw = 2.08 GeV) and the fact that /3 depends on the decay products 
of the decaying particle. For instance, the leptonic W decays involve only one charged lepton 
instead of two. As a result, we should use (|18D with [ln(M^/m^) — 1] for i = e,fi,T, 

which is scaled down by at least a factor of two compared with the Z-pair example. For 
W bosons decaying into an electron or positron, the resummed FSR distortion effects amount 
to a shift in the peak position of — 45MeV and a peak reduction factor of 0.86 per distorted 
resonance (i.e. 0.74 for a double-invariant-mass distribution), as can also be read off from (pOD.p| 

From the previous discussions it should be clear that FSR-induced distortion effects can be 
sizeable and should be taken into account properly in the Monte Carlo programs that are used 
for the ly-mass determination at LEP2. 
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